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ABSTRACT 

The influence of a linearized perturbation on stationary inflow solutions in an inviscid 
lO ! and thin accretion disc has been studied here, and it has been argued, that a pertur- 

bative technique would indicate that all possible classes of inflow solutions would be 
stable. The choice of the driving potential, Newtonian or pseudo-Newtonian, would 
1^ ' not particularly affect the arguments which establish the stability of solutions. It has 

then been surmised that in the matter of the selection of a particular solution, adop- 
lO ' tion of a non-perturbative technique, based on a more physical criterion, as in the case 

of the selection of the transonic solution in spherically symmetric accretion, would give 
' a more conclusive indication about the choice of a particular branch of the flow. 

I Key words: accretion, accretion discs - black hole physics - hydrodynamics - meth- 
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Oh' 
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^ ■ 1 INTRODUCTION 

The object of this work has been to study the stability of the stationary inflow solutions in axially symmetric accretion, under 
the influence of a linearized perturbation. The chosen model has been that of an inviscid and thin accretion disc. 
■ U se of this model, in which viscous effects have been ignored, has foun d a fairly widespread and consistent favour over the 

?H ' years ijAbramowicz fc ZtirekllT98ll : IChakrabartHllOSfll : iMolteni et al.lll99^ . As a consequence of such a physical prescription. 



there remains no mechanism for the outward transport of the angular momentum of the flow, and hence the specific angular 
momentum of the flow has been taken to be a constant of the motion. This naturally simplifies the analysis to a great extent, 
without significantly detracting from a true understanding of the underlying physi cs. The whole analysis has been based on 
the thin-disc configuration in accretion studies - a very well-known textbook model jPringlJlQSlllFrank et allll992l) that has 
been quite readily and regularly invoked by researchers in accretion astrophysics. 

It has been argued here, that in an inviscid and thin accretion disc, the stationary inflow solutions of abiding interest, 
are all stable under the influence of a linearized time-dependent pert urbation. The technique adopted has been sim ilar to the 
one already existing for the case of spherically symmetric accretion JPetterson et alJll980l : iTheuns fc Davidlll992ll . Separate 
emphasis has been laid on discussing the stability of the so-called transonic inflow solution - a solution that passes through 
the critical point (s) of the flow, and whose velocity is subsonic very far away from the accretor, becoming supersonic as the 
accretor is approached. The subsonic inflows, a class of physically meaningful inflow solutions, the maximum of whose velocity 
always remains below the velocity of the inflow solution passing through the critical point(s), have been quite extensively 
considered as regards their stability, since, like the transonic inflow, the subsonic infiow solutions also obey the same outer 
boundary condition that the flow velocity should decrease very far away from the accretor. And it may also be pointed out 
that since the subsonic solutions do not pass through any singular point of the flow, the mathematical problem of studying 
their stability is much simpler compared with that of the transonic solution. 

The classical Newtonian potential has been taken as the d riving factor which effects the infall of matter, although it 
is easy to see that the use of the pseudo-Newtonian potential of IPaczvnski fc Wiital il980l) - which satisfactorily represents 
general relativistic effects in a Newtonian framework and which is frequently in use to study the infall of matter under the 
gravitational influence of a black hole - alters nothing radically about the conclusions drawn from a linear stability analysis 
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of the inflow solutions. What it changes, however, is the number of critical points of the flow. In the context of the use of 
the Paczyriski-Wiita potential, mention may be made here of the fact that the number of critical points of the flow is also 
determined by various c onstants of the mo tion like accretion rate, specific angular momentum and the specific energy of the 
fiow. It has been shown l)Chakrabartilll99(t l that for a given angular momentum of the fiow, there exists a range of accretion 
rate, or equivalently a range of energy, such that the fiow has three critical points. Of these, the two outer points are saddle 
points, flanking between themselves, a third centre-type point. The flow, however, has only one critical point for values of 
accretion rate and energy which do not lie within the given range that will develop three critical points. 

I n solutions which admit of more than one saddle-type critical point, the existence of shocks is a possibility iChakrabarti 
I1990I) . For accretion on to a black hole, a solution would then pass through the outermost critical point (which is a saddle 
point), undergo a shock and then pass through the innermost critical point (again a saddle point) to cross the horizon with 
supersonic velocities. For accretion on to a neutron star, a flow passing through the inner critical point must undergo another 
shock to m eet the inner bound ary condition. The behaviour of these shocks are understood by invoking the Rankine-Hugoniot 
conditions (IChakrabartilll989l) . Once again it is to be emphasized that the outer boundary condition for all kinds of physically 
meaningful inflow solutions - shocked or continuous - is that the bulk velocity should decrease to zero very far away from the 
accretor. 



2 THE EQUATIONS OF THE FLOW 

For the radial drift velocity v and density p, the flow is governed by Euler's equation and the continuity equation. The radial 
component of Euler's equation is given by 

dv dv IdP dV „ 

where L, which is a constant, is the specific angular momentum, and the Newtonian gravity potential V = —GM/R. The 
choice of the pseudo-Newtonian potential of IPaczvhski fc Wiital l)l98ril would give V = —GM/{R — Vg), in which r-g is the 
Schwarzschild radius of a black hole. 

Use has been made of a polytropic equation of state 

P = kp'' (2) 

where 7 is the polytropic exponent, whose admissible range (1 < 7 < 5/3) is restricted by the isothermal limit and the 
adiabatic limit respectively. 

A s urface density E ha s been defined by vertically integrating p over the disc thickness H, and in the thin-disc approx- 
imation llFrank et al.l I ll992tl this gives T,=pH. The continuity equation, in terms of the surface density E, in the thin-disc 
approximation, is thus given by 

f + ij|M^» (3) 

Again invoking the thin-disc approximation, the disc thickness H is obtained from the vertical component of Euler's 
equation as 

Cs 

"(GM)i/2 

where Cs, the speed of sound, obeys the relation 

c,''=7V'- (5) 

The set of equations defines the whole problem completely. From Q and Q is obtained 

dv dv , dp dV 

m^^'dR^'^P M + dR-R^-' 

while lOl,© fi'iid © give 

I [p'^+^'/^j + R-'^'^ [p^''+'^/\R'^'] = (7) 
The two equations Q and Q above are those that govern the fiow in the inviscid and thin accretion disc. 



H^TT^R'^' (4) 



3 THE PERTURBATION EQUATION 

The steady-state solutions of (|SJ and ^ are vq and po, on which a re imposed small first-order perturbations v' and p' 
respectively. Following a similar prescription bv lPetterson et all (Il98d) for spherically symmetric accretion, it is convenient 
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to introduce a new variable, f=p^~'^^^^^vR^''^ , which in physical terms is the mass flux. Its steady-state value /o, as can be 
seen from (|7|, is a constant that is identified as the mass accretion rate. A linearized perturbation /', about /o is given by, 



f = Po 



{7-l)/2^5/2 



/ , 7+1 ' 

" Po H 7, — vop 



Linearizing in the perturbation variables v' and p', gives from l|7|, the expression. 



21 
dt 



-(7-l)/2^-5/2 



df 







7+1'" dR 
With the use of lO, successive differentiation of (jSJ with respect to time yields, first 



R' 



po 



dv' 
'dt 

and then 



-(7+i)/2 df ,df 



(7+l)/2 



dt 



dt^ ^''°dR[ dt 



Now linearizing in terms of the perturbation variables, gives from Q the equation, 
dv' d f , 2 p' \ ^ 



(8) 



(9) 



(10) 



(11) 



(12) 



where Cso is the speed of sound in the steady state. 

Differentiating 11211 with respect to time, and using the conditions given by |I5)|. II10II and lllll will finally deliver the 
perturbation equation as 

d\f' 



d 



dt^ 



dR 



df 



dt 



1 d 



vo dR 



Vo Vo 



2 2 

CsO 



7 + 1 



df 
dR 



= 



(13) 



This expression bears a very close resemblance in form to a perturbation equation obtained bv lTheuns fc DavidI 
using similar methods, for the case of spherically symmetric accretion. 



4 STABILITY ANALYSIS OF SOLUTIONS 

Before taking up the stability analysis of solutions, as implied by , it would be worthwhile to have an understanding of the 
critical points of the fiow. In the steady state, Q and Q can be combined to give the critical points - the required condition 
being that the numerator and the denominator of dvo/dR vanish simultaneously at those points. This leads to the expressions 



Vo = 



7 CsO 



7+1 
GM 



R^ 



7 + 1- R R^ ^''^ 

The second equation above, being a quadratic in R, indicates that there would be two critical points for the fiow. 
The choice of the pseudo-Newtonian potential of Paczvriski & Wiita (lOSOj), merely changes the number of critical points 
iChakrabartilll989t) . without altering the first of the two equations above. 

It is seen that the transonic infiow which originates very far away from the accretor, with a very low subsonic velocity, will 
pass through the first point - which is a saddle point, and then tend to curl around the second point - which is a centre-type 
point. In case the accretor is a large and distended object like an ordinary star (which justifies the use of the conventional 
Newtonian potential in this analysis), it is quite likely then that while the flow solution tends to curl around the inner critical 
point, it would also meet the surface of the accretor. This, as will be discussed later, has important implications for the 
stability of the transonic infiow. 

Returning now to 1131 1. a solution of the form /' — (;(i?)e~'"' is chosen, in which g(R) is the spatial part of the perturbation 
and Q is real. This leads to the result 



^9 ^, ^ d , , It/ 
-gn^-2in—{vog) + - — 
dR Vo dR 



vo[vo 



= 



(15) 



where =2/(7 + 1). 

For subsonic fiows, it is easy to understand that there are two values of i? - one close to the accretor and one very far 
away - where the perturbation, which is in the form of a standing wave, can be constrained to die out. Multiplying 1151 by 
vog and integrating over the range of R, which is bounded by the two points where the perturbation dies out, yields, 
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d7? 



(16) 



The above result could be obtained because the boundary terms all vanish. With a real value of g, for the subsonic flows - 
all of them governed everywhere by the condition vq'^ < P^Cao^ - fi^ will be positive, implying that the perturbation, which is 
in the nature of a standing wave, will be oscillatory an d at least will not grow in time. This line of reasoning was analogously 
established for spherically symmetric subsonic flows bv lPetterson et al] lIlQSd) . 

It is also possible to subject the subsonic flow solutions to an a nalysis where the pertu rbation may be treated as a travelling 
wave. The method for this, again follows a prescription used bv lPetterson et all lll98(J) for the spherically symmetric flow. 
The perturbation would have to be confined to a short-wavelength regime, in which the wavele ngth of the propagatin g waves 
will be much less compared to a characteristic length-scale of the system, which, in this case, as lPetterson et all lll980l) argued 
for the analogous spherically symmetric situation, could be the radius of the accretor, and that for an ordinary star is i?*. 
Since such a length-scale would be significantly large, it is somewhat justifiable to restrict the study to a perturbation of 
comparatively short wavelengths. Consequently, fi, the frequency of the waves, would be large. 

It is convenient to recast 1151 1 as 



dvo 2 1 d / 2\ r,.o 



dg 
dR 



dR 



5 = 



dR vo dR 

A solution to 11711 may be obtained by expanding g{R) as a power series JPetterson et al.|ll980l) . in the form 



g{R) = exp 



E 



k„iR) 



(17) 



(18) 



and to solve for the first few terms. 

For such a solution, terms in Q, and OP will respectively yield the conditions 

2 



dk-i 
dR 



dfc_i 
2iuo — TTT- -1 = 
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dR dR 
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dko 

Hr 



+ 



dvo 
dR 



f3^_d_ 
Vo dR ' 
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2ivo 
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(19) 



(20) 



(21) 



It is to be noted that as compared to the spherically symmetric case, the analogous results for the inviscid and thin 
accretion disc, given by 11911 . 1201 and 12H . vary by a scale factor of l3 for the speed of sound in the steady state. It is then 
quite straightforward to arrive at the expressions, 

^ = (22) 

dR vo±l3cso ^ ' 

and 

fco = — — In (/JuqCso) + constant (23) 

The first two terms in the power series expansion of g{R) are given by 12211 and 1231 1. As with the spherically symmetric 
case, here also self-consistency may be had by establishing 

n|fc-i(ii)| > > n-'|fci(i?)| (24) 

From the way it has been chosen, f2 ^ (i;o ± l3cso)/R*, it can be seen that the condition given by H24|l is satisfied, if the 
radial distance is not too small. For large R, the asymptotic behaviour of fci, fc_i and ko obtained from 121II . 1221 and 1231 1 
respectively, is given by 



k-i ~ _R , fco ~ In _R 



(25) 



The result given above upholds the self-consistency required by the condition in 1241 . The results in 1251 follow from a 
very possible power-law dependence of the steady flow velocity on the radius, at great distances from the accretor. 
The expression for the perturbation may then be written as 

sit 

' A+ exp ( ifl I — I -I- A- exp ( i^i 



/' 



^l3voCsO 



Vo + /3CsO 



dR 



vo - /3CsO 



(26) 
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In 1261 1 there is a linear superposition of two solutions with arbitrary constants A+ and A-, both travelling with velocity 
/Scso, relative to the fluid, which itself moves with a bulk velocity vo- One solution moves along with the bulk flow, while the 
other moves against it. 

Use of Q and ^ gives the perturbation in density as 



5/2^^-(^-l)/2 



V/3«oCi,o 



A, 



vo + (5CsO 



exp 



dR 



+ 



A. 



vo — (5CsO 



exp 



dR 



vo — (5CsO 



(27) 



vo + PCsO ^ 

Using 1271 1. the perturbation in velocity can then be obtained from and following some simple algebraic manipulation 
it can be rendered as 



/ , CsO p 
V =±— — 

P PO 

in which the positive and negative signs indicate outgoing and incoming waves, respectively. 
In a unit volume of the fluid, the kinetic energy content is 
1 



(po + p) {vo + v'Y 



(28) 



(29) 



The potential energy per unit volume of the fluid is the sum of the gravitational energy, the rotational energy and the 
internal energy, and is given by 

£'Pot^[po + p)—-[po + p)^+Poe + p-Q^^[pot)+-p g^(poe) (30) 

where e is the internal energy per unit mass jLandau fc Lifshitzlll987ll . 

The first-order terms vanish on time averaging. In that case, the contribution to the total energy in the perturbation 
comes from the second-order terms, which is given by 



tport = 2^0" +'"OpV+-p -Q^[PO^> 



(31) 



For an adiabatic perturbation, such as the one being prescribed here, the condition ds — can be imposed on the 
thermodynamic relation, de = Tds -\- {P/p^^dp, which will then give 



dpo' 



■(poe) 



CsO 

Po 



Combining the results from 1281 . 13H and 1321 will give 



po 



+ l) ± f3vo 



The expression for p' is to be obtained from 1271 . and upon time averaging, the result will be 



1 /3yl2 



Po 



-(7-l)/2 



(vo ± l3cso) 



CsO 

2 



+ l) ± Pvo 



(32) 



(33) 



(34) 



^ JO 

where fo = po^'^+^^^'^voR^^'^ , is a constant that in physical terms is the steady-state value of the accretion rate. The factor 
1/2 in 1341 ') derives from the averaging of p'^ . 

The total energy flux in the perturbation is obtained by multiplying fpert by the propagation velocity {vq ± l3cso) and 
then by integrating over the area of the cylindrical face of the accretion disc, which is 2n RH. Here H is to be substituted 
from ||1J and J^, and together with the fact that H ^ R in the thin-disc approximation iPrank et al.lll99'3l . an expression 
for the energy flux will be delivered as 

1/2 



fo \gm) 



± 1 



1-/3^ 



2P{M ± 13) 



(35) 



where M = vo/cso, is the Mach number. 

It is to be seen that for purely subsonic flows, M always remains less than which is a condition that is easily understood 
from the first of the two critical point conditions, given by 1141 . This would then imply that the total energy in the wave 
remains finite as the wave propagates. Indeed, for very great radial distances, when M — >■ 0, the flux is less compared to 
what it is near the outer critical point of the flow. The perturbation does not manifest itself as a runaway increase in the 
energy of the system. The system would then remain stable under the influence of such a perturbation. 

For the case of the transonic flow, a different line of reasoning would have to be adopted, since the condition vo"^ < /3^Cao^ 
would not hold everywhere, with the flow passing through a singular point. The physically meaningful flow in this case will 
be subsonic beyond the saddle-type outer critical point, pass through it to attain a supersonic velocity and then tend to curl 
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around the inner centre-type critical point, but not actually flow through it. Ixheuns fc DavidI if 992l) have established the 
stability of the transonic solution in spherical symmetry by arguing that the stability of the subsonic region of the flow, with 
the critical point as one boundary, will depend on conditions at the far outer end of the flow. In so far as the thin accretion 
disc is dependent only on R (the disc being invisci d and hence, the related equations are independent of any angle variable), 
the analogous argument of Ixhexms fc DavidI l)l992ll is expected to hold good to ensure the stability of the subsonic region of 
the transonic solution. 

For the supersonic region, which is finite, the reasoning of lGarlickl lIlQTSi) that a disturbance in this region will be carried 
away in a finite time, will also ensure the stability of the flow in this region. 

There is another important argument to uphold the stability of the transonic flow in a thin accretion disc. In the case of 
the spherically symmetric transonic flow, the velocity near the accretor has a n unbridled growth, su ch that any perturbation 
on it cannot be constrained to die out by a conceivable physical mechanism jPetterson et alJll980l) . In the inviscid and thin 
accretion disc, the transonic solution has no such runaway growth in the inner region of the flow, tending to curl around the 
inner critical point as it does. And while doing so, it may reach the surface of the accretor, if the accretor is as distended as an 
ordinary star. This can ensure that a perturbation on the solution in this region will not behave in an unconstrained manner. 

A note of caution has to be sounded here. The discussion presented above, pertai ns only to flows which a re continuous 
everywhere. However, for spherically symmetric flows with discontinuities like shocks, iTheuns fc DavidI JwO^) have argued 
that the expectation would also be one in favour of the stability of those flows. Under the present assumptions governing 
an inviscid and thin accretion disc, th e mathematical problem becomes very similar to the case of spherically symmetric 
accretion. In that case the argument of iTheuns fc DavidI il992l) for spherical symmetry could very likely be extended to the 
axially symmetric system being considered here. In any case, regardless of the geometry of the problem, it is a matter of 
general understanding that a disturbance in the supersonic region in a shocked solution, is transmitted completely through 
the shock and that in the subsonic region a disturbanc e is reflected at the sho ck, with a diminished reflected amplitude for the 
disturbance as compared with the incident am plitude (|Theuns fc Davidlll992l) . More to the point, by a local stability analysis, 
it has been demonstrated (|ChakrabartJll989ll that Rankine-Hugoniot shocks in an inviscid and thin accretion disc are stable 
under the influence of a short-wavelength perturbation in real time. 

It can then be safely established that a physically well-behaved and meaningful inflow solution is indeed stable under the 
influence of a linearized perturbation. This conclusion holds good for both the transonic inflow solution, as well as the entire 
class of subsonic inflow solutions. 



5 CONCLUDING REMARKS 

So far, it has been demonstrated that linear stability analysis, which is essentially a perturbative technique, offers no clue 
about the exclusive choice of a particular solution - indeed, to the extent that a perturbative technique can be relied upon, 
every solution (transonic or subsonic) seems just as much realizable as any other. It is then to be conjectured that the 
criterion for the selection of a particular solution would have to be non-perturbative in character and would have to be based 
on fundamentally physical arguments, like the maximization of the accretion rate or the minimization of the total energy 
associated w ith a solution, as in the case of the transonic inflow solution in spherically symmetric accretion jBondil 11953 : 
lGarlicklll97il . 

A recent work l|Rav fc Bhattacharieell2002l ) on spherically symmetric accretion has indicated that such an approach, 
however, would necessitate the addressing of the whole issue of the selection of a particular solution, from the viewpoint of 
the temporal evolution of the flow, instead of considering it solely on the basis of the stationary picture. The study has shown 
that it is the transonic solution which is decisively preferred to all the others. Such an expectation would also be valid for the 
case of the inviscid and thin accretion disc being discussed here. 

Support for such a contention can also be had from a very different quarter. Study carried out over the last two decades has 
established that there is a very close one-to-one correspondence between certain features of black hole physics and the physics 
of supersonic acoustic flows. More speciflcally, for an inviscid, barotropic and irrotational fluid flow (such as the conventional 
spherically symmetric accretion) the equation describing an acoustic disturbance can be rendered i n the form o f a metric 
that relates very closely to the Schwarzschild metric, represented in the Painlever-GuUstrand form (Ivisserlll997l ). Close to 
the sonic point of the flow, the similarity is particularly evident. Such a close correspondence between the unidirectionality of 
trajectories near a black hole and of flows passing through the sonic point, is a reason strong enough to suggest that transonic 
trajectories are the favoured ones. Extending that argument, a comparison may also be drawn for rotational, inviscid fluid 
flows, to have an understanding of the favoured status of the trajectories in this case. 
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